Auto-oscillation threshold and line narrowing in MgO-based spin-torque oscillators 
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We present an experimental study of the power spectrum of current-driven magnetization oscillations in MgO 
tunnel junctions under low bias. We find the existence of narrow spectral lines, down to 8 MHz in width at a 
frequency of 10.7 GHz, for small applied fields with clear evidence of an auto-oscillation threshold. Micromag- 
netics simulations indicate that the excited mode corresponds to an edge mode of the synthetic antiferromagnet. 

PACS numbers: 75.75.+a, 72.25.Pn, 85.75.-d 



Spin-transfer torques |1, 2], which involve the transfer of 
spin-angular momentum between spin-polarized currents and 
magnetization in magnetic multilayers, lead to many inter- 
esting possibilities for controlling magnetization dynamics 
in nanoscale devices. One important application concerns 
nanoscale radiofrequency (rf) oscillators, in which magneti- 
zation oscillations, driven by a DC spin-polarized current, are 
tunable with applied field and current with the potential for 
large power output. Such current-induced magnetization os- 
cillations have been observed in metallic spin-valve nanopil- 
lars Jsl] and nanocontacts 10] • However, there is only a limited 
body of work ||5l|^0,S] to date concerning clear spin-torque 
oscillations in magnetic tunnel junctions (MTJ) with large 
tunnelling magnetoresistance ratios (TMR), particularly for 
MgO-based junctions. An experimental realization of large- 
amplitude oscillations in such MTJ systems, which would al- 
low sufficiently large power to be generated for applications, 
is therefore an important scientific problem. 

In this article we present experimental evidence of magneti- 
zation auto-oscillations in MgO-based tunnel junctions which 
are in the "virgin" state, i.e. before noticeable deterioration 
of the tunneling characteristics due to large sustained applied 
voltages across the junction. We show that a clear auto- 
oscillation threshold is reached for increasing current, which 
is accompanied by a drastic increase of the emitted power 
and line narrowing. Furthermore, narrow spectral lines, with 
linewidths on the order of a few MHz for oscillation frequen- 
cies in the range of 8-15 GHz, are observed even at zero ap- 
plied fields. This is remarkable because high-frequency power 
spectra with comparable linewidths have only been observed 
in nanopillars until now in spin-valve systems under large ap- 
plied fields iSI^. As we will show, micromagnetics simula- 
tions indicate that the observed mode is an edge mode of the 
synthetic antiferromagnet. 

The measurements were performed on low resistance- 
area (RA) product MgO-based magnetic tunnel junc- 
tions with layer composition of Ta (3)/CuN (40)/Ta 
(5)/PtMn (20)/ CoFe3o(2) /Ru(0.8)/CoFe2oB2o(2)/Mg (1.3 
[nat ox])/CoFe2oB2o(3)/Ta(8) (thicknesses in nm) deposited 
on a Singulus Timaris tool. These junctions were patterned 



using electron-beam lithography and ion beam etching into 
rectangular nanopillars with a nominal size of 100x200 nm^. 
The electrical resistance of the tunnel junctions in the "vir- 
gin" state are typically 180 and 280 il for the parallel and 
antiparallel configurations, respectively, which correspond to 
a TMR ratio of 55%. The RA product is 0.9 n ^im^. Similar 
results were obtained with samples having an RA product of 
1.4 il/im^. The saturation magnetization hqMs is deduced 
to be 1.2 T in the patterned samples, while FMR measure- 
ments on unpatterned layer stacks show the effective magneti- 
zation to be 1 .9 T, the Lande factor to be 2. 14 and the uniaxial 
anisotropy to be 2.5 mT. Details concerning the characteriza- 
tion of the magnetic properties are presented elsewhere [10]. 

For the electrical measurements, the applied current 
through the tunnel junction is generated by applying a DC 
voltage to a 50 ohm resistor in series with the sample. Two 
bias tees allow DC and RF routing. One of the RF routes is 
amplified by 25 dB in a frequency band of 100 MHz — 26 
GHz and connected to a spectrum analyzer. The resolution 
bandwidth of the analyzer is set to 5 MHz for scans across the 
entire 0.1-26 GHz range, while zoomed scans are performed 
with a resolution bandwidth of 1 MHz. The final power spec- 
tra are obtained after subtracting a reference curve taken with 
a low current bias. For positive currents, electrons flow from 
the SAF system to the free layer, which favors a parallel mag- 
netization state in experiment. Material fatigue occurs for bi- 
ases above 0.4 V. In the study presented here, we use low bi- 
ases from —1.6 to +1.6 mA to study power spectrum in virgin 
state devices. The rf measurements are made in the presence 
of applied magnetic fields in the film plane, which are gener- 
ated with electrical coils or an electromagnet. In all our mea- 
surements a positive magnetic field along the easy axis favors 
the parallel magnetization configuration of the magnetic tun- 
nel junction. 

An example of a clear threshold behavior is presented in 
Fig. 1, where the power spectrum is shown for different ap- 
plied currents below and above the auto-oscillation threshold. 
The cuiTent threshold is estimated to be around / = 1.5 mA. 
When applying 7 = 1.5 mA through the device apeak appears 
at 10.59 GHz in the power spectrum that is well above the 
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FIG. 1 : Evidence of auto-oscillation threshold. Power spectral den- 
sity (PSD) of the voltage noise for applied currents from 1 .4 mA to 
1.6 mA. Inset: Power spectrum for / = 1.6 mA over the entire mea- 
surement bandwidth of 0.1-26 GHz. 



noise floor, in contrast to a relatively flat spectrum for / < 1 .5 
mA. If the current is increased further above the threshold 
value, a large increase in the total power is measured and the 
linewidth at half maximum reduces to a few MHz, with 8 MHz 
being the narrowest linewidth observed at room temperature 
and under the experimental conditions described. These spec- 
tra were measured at zero applied magnetic fields but after 
preparing the magnetic tunnel junction in the antiparallel (AP) 
configuration, i.e. between every measurement we applied a 
strong field to saturate the free layer magnetization along the 
negative easy-axis direction. The frequency of the peak is not 
observed to vary much with the current above threshold, in 
contrast to free-layer excitations above threshold for which a 
large redshift is expected for in-plane fields. The influence of 
the current can also be seen in the Stoner-Wohlfarth astroid 
for in-plane fields as a strong distortion occurs under positive 
cuiTent with respect to a low current astroid. The positive cur- 
rent leads to an additional torque that favors the parallel (P) 
state. At very low frequencies, a large contribution from 1// 
noise is observed. At low current thermal modes are also vis- 
ible (not shown). From the variation of the frequency of the 
thermal FMR mode as a function of applied field, we deter- 
mined the effective saturation magnetization Mcs from fits to 
the Kittel formula, uj = iy/{Ha + Hk){Ha + Hu + Mcs), 
where Ha is the applied field and Hk is the anisotropy field, 
from which we find Afcff = 1.2 T. 

The full measured power spectra across the range of 0. 1-16 
GHz, as a function of applied field along the easy and hard 
axes, are presented in Fig. |2]as a color map. In scans along 
both field directions, we observe a number of peaks in the 
power spectrum but only one peak attains large power (> 100 
nV^/Hz), which corresponds to the peak shown in Fig. [T] A 
second harmonic of this peak is also observed (as shown in the 
inset of Fig. [U, around 20-22 GHz, but this is not presented 
in the color map for the sake of clarity. Third harmonics were 
not seen as the frequencies are beyond the scope of the ap- 
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FIG. 2: (Color online) Density plot of the power spectrum as a func- 
tion of applied magnetic field along the (a) easy axis and (b) hard 
axis, for an applied current of / = 1.6 mA. 



paratus. It is interesting to note that the main peak persists 
and retains a large power even at zero applied field, which 
can easily be seen for field sweeps along both the easy and 
hard axis directions. Moreover, the intensity of the peak is ob- 
served to increase as the applied field magnitude is decreased. 
Other peaks corresponding to thermally-activated modes, in 
the range of 5-7 GHz, are present over a wide range of the 
fields studied. A more detailled discussion of these modes 
will be presented elsewhere. 

Micromagnetic simulations show that the observed modes 
coiTespond to edge modes of the synthetic antiferromagnet 
bilayer. The simulations consist of tracking the time evolu- 
tion of magnetization of the full free layer/synthetic antifeiTO- 
magnet/antiferromagnet system using the same materials, film 
thicknesses and TMR ratio as the experimental system. In our 
notation for the synthetic antiferromagnet, SAFl refers to the 
magnetic layer adjacent to the free layer and SAF2 refers to 
the magnetic layer pinned by the PtMn antiferromagnet. The 
system is discretized with a finite element method with a lin- 
ear basis function and a discretization size of 2 nm, which is 
below the exchange length of 4.5 nm flT]. Three different 
calculations were performed for a given applied cuiTent and 
applied field geometry. In the first, a mutual spin-torque is 
applied to all magnetic layers self-consistently. In the second, 
the spin-torque is applied only on the free layer, and in the 
third, the spin-torque is applied only on SAFl. This allows us 
to clearly identify the excitations in the spin-torque-coupled 
free layer/SAFl layer system, the free layer, and the SAFl 
layer separately. The threshold current for oscillations in the 
simulations, 6 mA, is higher than the observed experimen- 
tal value of 1 .5 mA, which can be explained by differences in 
spin polarization or other transport parameters. For an applied 
field of 5 mT along the easy direction we observe SAFl layer 
excitations at 11.5 GHz and 23.1 GHz (second harmonic), as 
shown in Figure [3] In a second simulation run we applied an 
external field of 5 mT in the hard axis direction and observed 
a redshift in the frequency, in contrast to the behavior for easy 
axis field variations and in agreement with the experimental 
observations. For applied fields along both the easy and hard 
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FIG. 3: Simulated power spectral density of edge modes in the syn- 
thetic antiferromagnet, for an applied field of 5 mT along the (a) 
easy and (c) hard axes. The micromagnetic configuration of the top 
synthetic antiferromagnet layer (SAFl), illustrating the excited edge 
mode, are shown in (b) and (d) for the easy- and hard-axis fields, 
respectively. Good quantitative agreement is found between the sim- 
ulated and experimental frequencies of the main excitation peak. 
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FIG. 4: (a) Frequency, (b) linewidth (squares) and integrated power 
(diamonds) of the main peak as a function of applied field along the 
easy axis, with / = 1.6 mA. 



axis directions we find good quantitative agreement with the 
experimental frequencies. 

The variation of the frequency, linewidth, and total power of 
the experimental peak, as a function of the easy axis field un- 
der an applied current of 1 .6 mA, is presented in Fig.lD These 
parameters were extracted using a Lorentzian fit to the power 
spectra. As we mentioned above, the linewidth decreases as 
the applied field is reduced. Peaks observed in this low-field 
region are much narrower than those at higher fields, with 
a minimum in the linewidth of 8 MHz observed at zero ap- 
plied field, after which the line broadens again as the switch- 
ing transition is approached. Such narrow lines are consistent 
with the auto-oscillation regime IT3I1 . Some jitter in the 
spectral lines was observed by measuring the spectra at differ- 
ent sweep frequencies, which indicates that inhomogeneous 
broadening may be an important factor for narrow lines and 
that the actual oscillator coherence is much better than what 
the measured linewidths suggest. Further experimental stud- 



ies on this issue are currently in progress. The maximum total 
power is observed at low applied field, which indicates that the 
oscillation, in addition to better coherence, also gains in am- 
plitude. A similar behavior for the linewidth and power was 
observed for field sweeps along the hard axis, with the excep- 
tion of a much smaller window of existence for the excited 
mode, which can easily be ascertained from Fig.|2]3. 

In summary, we provide clear experimental evidence of 
current-driven magnetization oscillations above the auto- 
oscillation threshold in magnetic tunnel junctions based on 
MgO barriers. We have brought to light the possibility of ob- 
taining narrow spectral lines, down to 8 MHz at a frequency 
of 10.7 GHz, at zero applied magnetic fields, in tunnel junc- 
tions close to their "virgin" state. Micromagnetics simulations 
show that the excited mode is the edge mode of the synthetic 
antiferromagnet bilayer. 
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